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IN TRODU C TION.  The  main  concern  of  the  Land  Locomotion  Lab- 
oratory, ATAC,  is  the  relationship  of  a vehicle  to  the  terrain  over  which 
it  travels.  Once  an  insight  into  this  relationship  has  been  attained,  the 
problems  encountered  from  the  initial  design  of  a vehicle  to  its  ultimate 
use  in  the  field  can  be  more  easily  grasped  and  rationally  solved.  The 
solution  of  engineering  problems  depends  upon  the  selection  of  the 
relevant  variables  and  a description  of  the  functional  relationship  among 
these  variables.  The  selection  of  the  vehicle  variables  are  within  broad 
limits  at  the  disposal  of  the  designer,  but  for  the  terrain  or  soil,  the 
selection  of  suitable  variables  becomes  much  more  complicated.  Soil 
is  probably  one  of  the  most  complex  of  all  engineering  materials  £lj  . 
Researchers  in  soil  mechanics  have  added  much  to  the  knowledge  of  the 
mechanical  and  physical  characteristics  of  soils,  but  as  yet,  no  fully 
satisfactory  general  theory  is  available. 

Land  locomotion  is  an  engineering  application  of  soil  mechanics  to 
off-road  vehicular  operation.  Its  objective  is  to  determine  the  relation- 
ships of  vehicles,  or  more  precisely,  of  the  wheel  and  track,  to  the 
strength  properties  of  the  soil.  In  land  locomotion  research,  one  of  the 
important  questions  is  the  nature  of  the  pres  sure -sinkage  relationship. 
The  Bekker  equation 


represents  a family  of  curves  with  three  unknown  constants,  kc,  k^, 

and  n,  and  two  variables,  p the  pressure  under  a loaded  area,  and  z 
the  sinkage*  J^2j  . 

* b is  a known  constant,  the  plate  width.  p is  per  unit  area  (i.  e.  , 
square  inch). 
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A set  of  these  constants  which  will  approximately  describe  pres sure  - 
sinkage  observations  can  be  obtained  Q3j  . To  determine  these  constants, 
pressure-sinkage  experiments  were  performed  in  the  laboratory  with 
footings  of  various  sizes  and  shapes.  The  constants  obtained  were  used 
to  predict  the  sinkages  for  other  loaded  areas.  It  was  found  that  the  pre- 
dictions were  adequate  for  tracked  vehicles,  where  the  relative  shapes 
of  the  loaded  areas  were  similar.  Rectangular  test  footings  with  a 
length/width  ratio  greater  than  5 were  used  for  determining  the  soil 
parameters  used  in  the  prediction  equations.  The  basic  equation  in- 
cludes only  the  width  term  b.  For  cohesionless  soils,  the  ultimate 
bearing  strength  is  dependent  on  the  width  only  for  long  loaded  areas  (.4) . 
To  determine  what  the  minimum  length  had  to  be  before  a footing  was 
not  considered  long,  laboratory  tests  were  conducted  with  footings  of 
varying  ~tf b ratios.  Acceptable  results  for  the  pressure-sinkage 
relationship  were  obtained  when  the  -2/b  ratio  was  greater  than  5. 
Consequently,  all  pressure-sinkage  measurements  were  taken  with 
footings  of  at  least  an  X/b  ratio  of  5 or  greater. 

The  ultimate  equations  in  which  these  soil  strength  parameters  are 
to  be  used  apply  to  the  general  case  of  predicting  vehicle  performance 
for  both  tracked  and  wheeled  vehicles.  As  noted  above,  the  predictions 
of  tracked  vehicle  sinkage  and  motion  resistance  have  been  generally 
satisfactory.  For  wheeled  vehicles,  however,  improvements  are 
needed.  One  of  the  differences  to  be  noted  between  a tracked  and  wheeled 
vehicle  is  the  shape  of  the  loaded  area.  In  most  caBeB,  tracked  vehicles 
have  a contact  area  of  relatively  long  length  as  compared  to  width. 

Such  a length-width  ratio  i6  not  the  usual  situation  for  wheeled  vehicles 
at  moderate  sinkages.  Most  tires  will  have  a contact  area  where  the  4 /b 
ratio  is  close  to  1 or  2.  * Therefore,  it  was  thought  that  the  shape 
of  the  loaded  area  when  X/b  was  less  than  5 might  have  effects  on 
the  pressure-sinkage  relationship.  A clearer  understanding  of  the 
pressure-sinkage  relationship  in  this  region  would  provide  us  with  an 
improved  soil-vehicle  model  with  broader  and  more  useful  applications. 
Consequently,  a test  program  was  undertaken  by  the  Land  Locomotion 
Laboratory  to  investigate  further  the  pressure-sinkage  relation. 


#For  example,  the  Army  6x6  5-ton  truck  normally  carries  an 
11.  00  x 20  tire.  At  one  inch  sinkage,  I b - 1 and  at  eight  inch 

sinkage  Alb  = 2.33. 
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DESCRIPTION  OF  TEST  PROGRAM.  The  test  program  to  study 
the  effect  of  plate  size  on  the  vertical  soil  strength  or  sinkage  para- 
meters was  divided  into  two  parts.  The  first  part  comprised  a study 
of  the  reliability  or  repeatability  of  the  test  results  for  the  equipment  and 
mixing  techniques  that  were  to  be  used  in  the  tests.  The  second  part 
covered  the  measurement  of  the  load- sinkage  curves  while  using  differ- 
ent sized  plates. 


In  carrying  out  the  experiments  a laboratory  model  bevameter  is 
used  to  drive  a constant  speed  hydraulic  piston  arrangement  which  pushes 
a plate  into  the  soil  in  a bin.  The  depth  of  sinkage  and  the  force  on  a 
load  cell  are  plotted  electrically  by  an  X-Y  plotting  device.  The 
curve  is  traced  on  linear  graph  paper.  Prom  this  graph  values  for  p 
and  z may  be  read  off  and  plotted  on  log-log  paper.  The  slope  of  the 
least  square  fitted  line  on  the  double -log  plot  gives  an  estimate  of  the 
parameter  n.  The  constant  term  in  such  a fitted  equation  is  the  logarithm 
of 


where  b is  the  width  of  the  plate  used.  Thus  it  is  seen  that  kc  and  k^ 
cannot  be  estimated  by  any  straight  forward  statistical  technique.  Use  of 
two  different  b values  .however,  will  permit  settingup  two  simultaneous 
equations  in  kc  and  k^ . 

PART  I.  We  wished  to  determine  the  maximum  number  of  pene- 
trations that  could  be  made  with  one  preparation  of  the  soil  bin.  For 
this  purpose  we  set  up  a uniformity  trial  using  only  a 2"  x 10"  plate. 
Orientation  and  location  of  the  penetrations  was  arranged  as  indicated 
in  Figure  1.  The  three  orientations  shown  were  carried  out  in  a ran- 
domized block  design  with  six  replicates. 
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j 


FIGURE  I:  SAMPLE  BLOCK  LAYOUT  FOR  UNIFORMITY  TRIALS 

RESULTS  FOR  PART  I:  Mean  pressures  were  computed  for  the 

\ orientations  at  1",  2"  and  3"  depths  of  sinkage.  Mean  differences  for 

* 

the  orientations  at  a given  depth  were  found  to  be  homogeneous  (i.  e.  , 
not  significant).  At  the  2"  depth  of  sinkage,  the  coefficient  of  vari- 
ation was  about  6 percent,  quite  a sati  s factory  value.  Examination 
of  the  variability  within  the  orientations  showed  that  the  variation  within 
the  arrangement  A^  was  significantly  greater  than  that  for  the  other 
arrangements. 

CONCLUSIONS  FOR  PART  I:  The  experimental  procedure  would 
yield  results  of  adequate  reliability.  The  orientation  Aj  appeared  un- 
desirable for  taking  pressure-sinkage  measurements.  Therefore,  we 
decided  to  make  the  spacings  between  determinations  as  similar  to  the  A 
arrangement  as  practicable.  Since  a rectangular  plate  of  size  3”  x 10" 
was  to  be  used  in  Part  2,  it  appeared  necessary  to  "beef  up  the  apparatus" 
to  handle  the  greater  pressures  required  to  sink  such  a large  plate. 

PART  2.  The  second  phase  of  the  test  program  comprised  measure- 
ment of  the  load- sinkage  relation  in  dry  sand  with  plates  of  varying  sizes. 
Analysis  of  test  data  would  provide  estimates  of  the  three  parameters  in 
j the  Bekker  equation  and  their  variations  for  the  sizes  and  shapes  of 

j plates  used.  The  results  should  show  the  dependence,  if  any,  of  the 

J parameters  on  the  size  and  shape  of  plate. 

Many  problems  arose  in  the  consideration  of  the  second  part  of  the 
test  program.  One  point  was  that  all  plates  should  be  tried  in  one  mix 
of  the  bin.  Thus,  a complete  block  design  would  be  preferred. 
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A second  point  is  that  the  sinkage  equation  is  a stress-strain 
relationship.  It  may  also  be  described  as  a functional  relation  Cb]  . 

The  problems  of  estimation  which  arise  for  the  functional  relation  have 
been  resolved  by  Dr.  Joseph  Berkson  by  introducing  the  ’’control  vari- 
able" concept  . It  appears  that  z may  be  taken  as  the  control  vari- 
able in  our  problem.  This  approach  is  contrary  to  the  usual  dependent- 
independent  variable  point  of  view,  but  Berkson  has  shown  that  the  method 
is  unbiased  for  estimating  the  functional  relation  if  the  errors  in  p are 
unbiased. 


Thirdly,  it  was  clear  that  a statistical  analysis  of  the  estimation  pro- 
cedure for  the  Bekker  equation  was  needed.  With  transformation  to  loga- 
rithms of  p and  z it  is  assumed  that  the  standard  linear  regression 
assumptions  are  valid  in  the  transform  £7}  • Thus,  estimation  of  the 
parameter  n is  quite  straightforward.  When  n has  been  obtained, 
the  procedure  takes  z = 1,  hence,  log  z = 0,  and  predicts  a value  of 
log  p,  say  Pq.  Now,  anti-log  PQ  = p#  = k ^ + kc/b.  By  taking  two 
values  of  the  plate  width,  bj  and  b^,  and  corresponding  and  p^* 

values,  the  estimation  equations  for  kr  and  k^  become 


(2) 


k = 
c 


(Pl*  - P2*)bjb2 


b2  - bl 


and 

(3) 


<P 


b2p2*  " biPi* 


These  results  pointed  out  two  things: 

1.  Widely  spaced  b values  to  permit  use  of  large  b -b. 
reduce  the  variance  of  the  estimates  of  kQ  and  k^  , 

2.  A formula  for  the  variance  of  p#  is  needed. 


would 


The  variance  formula  for 
a nonlinear  function  of  P, 


p#  presents  some  difficulty  because  p*  is 
. We  recall,  however,  that  the  choice  of  an 


appropriate  experimental  design  will  give  us  direct  estimates  of  the 
experimental  error  for  our  estimated  k_  and  ki  values  so  that  we 


can  bypass  the  vari  ance  formula  problem. 


<P 
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DESIGN  OF  THE  PART  2 EXPERIMENTS.  From  the  statistical 
analysis  it  was  clear  that  the  largest  possible  difference  in  plate  widths 


should  be  used  to  estimate  the  parameters  kc 


This  consider- 


ation along  with  the  desire  for  a complete  block  experiment  already  men- 
tioned resulted  in  a revision  of  the  choice  of  dimensions  for  the  rectangu- 
lar plates  to  be  used  in  the  experiment.  The  plate  sizes  actually  selected 
are  given  in  Table  1. 


TABLE  1 

LIST  OF  PLATE  DIMENSIONS  FOR  PART  2 EXPERIMENTS 


(in  inches) 


Rectangles: 


Circles: 

Diameters: 


1x4 

1x6 

1x8 

1 x 10 

2 and  4. 


2x4 
2x6 
2x8 
2 x 10 


3x4 
3x6 
3x8 
3 x 10 


The  entire  set  of  plate  s thus  provided  14  treatments  to  be  setup 
in  a completely  randomized  block  experiment  where  one  block  equals  a 
mix  of  the  soil  bin.  Each  plate  was  to  be  used  in  a randomly  selected 
plot  of  size  about  2'  x 2'  to  make  a single  measurement  of  the  pressure' 
sinkage  curve.  It  was  further  decided  to  complete  six  replicates  which 
would  yield  24  independent  pairs  of  estimates  of  kc  and  k^  , but 
would,  of  course,  yield  84  estimates  of  the  parameter  n,  one  estimate 
being  obtained  from  each  pressure  - sinkage  curve. 

ANALYSIS  OF  RESULTS  FOR  PART  2.  Analysis  of  the  estimated 
values  for  n is  straightforward  and  results  are  readily  interpreted. 

For  the  kc  and  k^  values  some  difficulties  arise.  Plotting  the  means 
is  a most  useful  device  for  aid  in  understanding  the  effects  indicated  by 
the  analysis. 

Figures  2 and  3 show  the  width  and  length  effects  on  n without 
considering  the  interaction.  In  order  to  present  the  interaction  effect, 
we  show  the  usual  type  of  plot  for  displaying  an  interaction.  Parallel 
lines  for  the  various  lengths  of  plate  would  be  indicative  of  no  interaction. 
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Thus,  the  lack  of  parallelism  exhibited  in  Figure  4 is  indicative  of  the 
nature  and  source  of  the  interaction.  The  major  pattern,  however,  is 
still  that  shown  in  Figures  2 and  3.  There  appears  to  be  a decrease  in 
the  n values  with  an  increase  in  area  of  plate  whether  the  area  increase 
is  due  to  change  in  length  or  change  in  width.  The  width  effect  is  much 
smaller  in  going  from  2"  to  3",  but  the  length  effect  is  about  the  same 
at  all  widths. 

In  considering  the  analysis  of  the  estimated  kc  and  k<£,  values  it 

will  be  useful  to  recall  equations  (2)  and  (3)  which  show  how  the  estimates 
are  obtained.  There  are  some  obvious  points  to  be  noted  from  these 
equations,  but  we  shall  defer  them  until  later.  Each  experimental  unit 
or  plot  in  the  soil  bin  yields  a p versus  z relation  as  drawn  by  the  x,  y 
plotter  while  a single  plate  is  sunk  into  the  soil.  Then  results  from  two 
different  plate  widths  have  to  be  combined  to  obtain  a single  estimate  of 
kc  or  k <4>  . We  may  combine  1"  and  2”  widths  or  widths  of  2"  and  3" 
or  1"  and  3".  As  shown  above,  the  latter  is  the  best  choice,  but  we 
see  that  a single  replicate  or  "set"  will  only  give  us  four  such  indepen- 
dent estimates. 

By  keeping  the  width  difference  constant  and  varying  the  length  we 
can  pair  within  one  set  these  four  pairs  of  plates: 


1x4  and  3x4  1x8  and  3x8 

1x6  and  3x6  1 x 10  and  3 x 10 


to  give  us  four  pairs  of  estimates  for  kc  and  k . Utilizing  the  six 
"sets"  gave  us  24  values  for  analysis.  * Estimates  obtained  from  these 
pairings  may  be  analyzed  for  the  length  effect  alone.  But  we  would  also 
like  to  study  the  width  of  plate  effect  on  these  two  parameters  and  look 
for  interaction,  if  any,  as  we  did. for  the  parameter  n. 

First,  we  consider  the  estimates  of  kc.  The  mean  values  obtained 
for  the  pairs  just  listed  were: 


& 

A set  consists  of  six  replicates  carried  out  at  the  same  depth.  Eight- 
teen  replicates  were  actually  completed  at  each  of  two  depths.  Aggre- 
gation over  six  replicates  forms  a set.  Thus,  there  are  a total  of  six 
sets. 


ESTIMATED  n VALUES,  BEKKER  PRESSURE-SINKAGE  EQUATION 
rectangular  plate,  dry  sand 
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Figure  2 
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ESTIMATED  n VALUES  , BEKKER  PRESSURE  - SINKAGE  EQUATION 
rectangular  plates,  dry  sand 


WIDTH  OF  PLATE  (INCHES) 


WIDTH  OF  PLATE  (INCHES) 
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Pairs: 


1 x 4 
and 
3x4 


1x6 

and 

3x6 


1x8 

and 

3x8 


1 x 10 
and 
3 x 10 


-3.  305 


Mean  Value  of  kc 
-4.030  -4.000  -4.166 


Relevant  mean  squares  abstracted  from  the  analysis  of  variance  C 8 J for 
comparing  these  means  are: 


Source 

d.  f. 

M.  S. 

Length  of  Plate 

3 

0. 8977 

Length  of  Plate  by  Depth 

3 

1.  0526 

Length  of  Plate  by  Sets 
within  Depth 

12 

0.  2262 

Since  the  P(F  3.  89)  = 0.  05,  we  may  conclude  that  there  is  a 

J I lb  ,T 

length  effect  on  kc  and  that  the  length  effect  is  not  uniform  at  the  two 
depths  studied  (i.  e.  , F = 3.  97  = 0.  89770.  2262). 

To  extend  our  analyses  other  plate  pairings  were  studied.  Finally, 
it  seemed  in  our  judgment  that  the  following  pairs  were  most  useful  for 
securing  information  on  the  width  effect  plus  a little  more  length  infor- 
mation: 


2x4 

and 

3x8 


2x6 

and 

3 x 10 


1 x 4 
and 
2x8 


1x6 

and 

2 x 10 


Pairs: 
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Mean  Value  of  kc 


-2.  463  -2.495  -5.  169  -5.063 


Note  that  the  pairs  selected  show  a constant  length  difference  of  4 inches. 
Furthermore,  the  length  differences  are  balanced  for  comparing  1 and 
2 inch  widths  with  2 and  3 inch  widths  for  the  study  of  the  plate  width 
effect  in  estimating  kc>  It  appears  that  the  latter  effect  is  large. 

Relevant  mean  squares  from  the  analysis  of  variance  are  as  follows: 


Source 

d.  f. 

M.  S. 

Width  Pairs 

1 

41.  717 

Width  Pairs  by  Depth 

1 

7.  183 

Length  Differences  within  Width  Pairs 


Length  in  1 and  2 

1 

0.003 

Length  in  2 and  3 

1 

0.034 

Length  Difference  by  Depth 

in  1 and  2 

1 

0.  006 

in  2 and  3 

1 

5.  914 

Experimental  Error 

12 

1. 529 

For  comparing  the  width  pairs  we  obtain  an  F ratio  = 41.717/1.  529>27. 

We  conclude  that  the  parameter  kc  depends  on  the  choice  of  plate  widths. 
In  estimating  kc,  pairing  of  1 and  2 inch  widths  gives  a larger  value 
(algebraically)  than  pairing  2 and  3 inch  widths.  Looking  at  it  another 
way  the  data  show  that  changing  the  length/ width  ratio  for  the  plates  from 
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the  range  4 to  6 to  a range  of  2 to  3 materially  alters  the  kc  value 
obt  ained. 


Similar  analyses  for  the  estimates  of  the  parameter  k were  carried 
out  as  just  described  for  kc.  We  note  that  the  k^>  estimates  were  posi- 
tive, mostly  in  the  range  +7  to  +10.  Our  conclusions  about  the  effect 
of  plate  width  on  the  parameter  k were  the  same  as  reached  in  regard 
to  kc,  but  no  length  effect  was  detected. 

With  these  results  available  we  may  return  to  the  further  consideration 
of  the  method  used  to  obtain  the  estimates.  First,  we  observe  that  if  there 
were  no  width  of  plate  effect  in  the  p versus  z relation,  then  k<jb  would 
be  equal  to,  say,  p*.  In  the  notation  used  above,  no  width  of  plate  effect 
would  mean  P01  = P03’  if  one  inch  and  three  inch  plate  widths  were  used, 
and,  hence,  p-j#  = p^*  = a common  value,  p*.  Observed  p-j*  and  p^* 

values  would  differ,  of  course,  due  to  experimental  variation,  but  average 
value s would  be  equal.  Further,  when  k^j  = p*,  then  the  solution  for  kc 

is  zero.  The  parameter  kc  was  once  regarded  as  a measure  of  cohesive - 
ness  of  the  soil  and,  hence,  might  be  zero  for  sand  [3j  . The  present 
experiments,  in  which  mason's  sand  was  used  certainly  do  not  agree  with 
this  point  of  view  about  the  parameter  kc> 


Recognizing  then  that  the  present  model  for  the  p versus  z relation 
does  admit  a width  of  plate  effect  our  analyses  so  far  have  indicated  the 
magnitude  of  this  effect.  In  addition,  we  have  obtained  some  indication  of 
a length  effect  on  kc  although  none  appeared  for  k-L  . 


Second,  with  respect  to  the  method  of  obtaining  the  estimates  of  kc 
and  k <p,  we  point  out  that  the  estimates  obtained  are  correlated.  This 
correlation  cannot  be  avoided  because  the  values  are  based  on  the  solution 
of  two  simultaneous  equations.  In  fact,  the  estimates  also  would  have  been 
correlated  if  we  had  been  able  to  obtain  them  by  a direct  least  squares 
procedure.  In  the  least  squares  case,  however,  it  is  usually  easy  to  write 
down  the  covariance  of  the  estimates  and,  hence,  to  find  the  correlation 
if  desired.  For  the  two  simultaneous  equations  in  k£  and  k^  , matrix 

methods  can  be  applied  to  find  the  covariance  matrix  for  kc  and  k Q . 

The  re  suit  ob  tained  for  the  covariance  of  kc  and  k^is 
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(4) 


-1 


cov(kc,  k<£)  = (b  _b  )2  (bj  b2  V<PiS‘>  +bib2  V(P  *)) 

^ I “ 


where  b£  is  the  width  of  the  larger  plate  and  V(  ) denotes  variance  of 
the  enclosed  quantity. 

It  is  interesting  to  examine  the  consequences  for  different  choices  of 
plate  widths  in  estimating  kc  and  k^.  Asa  first  approximation  it  may 

be  adequate  to  assume  V(p^*)  — V(p2*).  If  this  assumption  is  made, 
the  following  correlation  values  are  obtained: 


Plate  widths  paired:  (widths  in  inches) 

1 and  3 1 and  2 

Value  of  correlation  of  k 

c 

-0.8944  -0.9487 

Observed  values  of  the  correlation  between  kc  and  k0  will,  of  course, 

differ  from  these  theoretical  values  because  of  experimental  variation  and 
lack  of  equivalence  of  V(p^*)  and  V(p£*).  It  is  clear,  however,  that  the 

above  results  further  support  the  use  of  maximum  difference  in  plate  widths 
for  estimation  of  these  parameters.  Actual  plots  of  the  k£  and  k pairs 

support  the  high  correlation  of  these  estimates. 

These  considerations  and  the  analyses  for  k£  and  k^p  made  it 
seem  desirable  to  return  to  an  analysis  of  the  84  p*  values  which  corres- 
pond to  the  84  n values  analyzed  earlier.  Such  an  analysis  should  give  a 
clearer  assessment  of  the  length,  width,  and  interaction  effects  than  we 
have  obtained  from  analysis  of  the  kc  and  k^>  values.  Any  effects  found 

in  this  analysis  must  then  be  present  in  the  k£  and  k ^ values  because 
of  the  method  of  derivation  of  these  estimates  from  the  p*  values. 


2 and  3 


and  k 


4> 

-0.  9805 
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The  means  for  the  p#  values  are  presented  in  Table  2.  Table  3 
following  gives  the  analysis  of  variance  for  the  p*  values.  From  the 
table  of  means  it  is  seen  that  the  width -of-plate  effect  is  much  larger 
than  the  length -of-plate  effect.  Both  effects,  however,  may  be  judged 
significant  from  the  analysis  of  variance  results.  It  is  interesting  to 
note  that  the  length  by  width  interaction  mean  square  is  so  small;  no  in- 
teraction is  indicated.  This  result  is  in  contrast  to  the  analysis  of  the  n 
values,  Table  4,  for  which  the  interaction  effect  was  judged  significant. 

CONCLUSIONS  FOR  PART  2.  A large  series  of  experiments  have 
been  conducted  to  study  the  pressure  versus  sinkage  relation  and  to  esti- 
mate the  parameters  in  this  relation.  The  experiments  were  carried  out 
using  a relatively  homogeneous  soil  material,  dry  mason's  Band.  Bearing 
plates  used  comprised  two  circles  of  2 and  4 inch  diameters  and  12 
rectangular  plates  varying  in  size  from  I x 4 to  3 x 10  inches 
(refer  Table  1 for  list  of  plate  sizes.  ) Thus,  the  length  over  width  ratio 
of  the  plates  ranged  from  a maximum  of  10  to  1 down  to  4 to  3,  or 
from  long  narrow  plates  to  almost  square  plates.  This  range  of  length 
over  width  ratios  was  selected  to  cover  the  range  from  tracked  vehicles 
to  wheeled  vehicles  with  tires. 

From  the  estimates  of  kc>  k^  , and  n plus  the  estimated  pressures 
for  one  inch  sinkage  our  analyses  show  that: 

( 1)  The  parameter  n varies  with  the  width  and’ length  of  plate.  For 

- - one  inch  width  of  plate  the  n value  was  0.  83;  at  three  inch  width, 

the  value  was  0.  68.  Between  2"  and  3”  width  there  was  little 
change.  With  length,  n varied  from  0.78  to  0.70  over  lengths 
of  4"  and  10".  The  decrease  was  nearly  uniform  over  the  6" 
inte  r val. 

(2)  The  estimated  pressures  for  one  inch  of  sinkage  (the  p*  values) 
show  variation  with  both  length  and  width  of  plate  but  no  interaction 
of  the  factors  is  indicated.  (Refer  to  Tables  2 and  3.  ) 

(3)  The  parameter  kc  decreases  algebraically  with  increase  in  plate 
length.  The  algebraic  change,  however,  was  much  greater  when 
estimates  were  compared  from  pairing  of  plates  of  1"  and  2" 
widths  with  estimates  obtained  by  pairing  2"  and  3"  widths. 

(4)  The  parameter  k^>  showed  little  or  no  response  to  length  of 
plate  but  a large  response  to  width  of  plate. 
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(5)  The  estimates  of  kc  and  are  highly  correlated.  This 

correlation  is  negative  so  that  when  kc  increases  in  value  k^ 
decreases  in  value. 


From  these  analyses  it  appears  that  the  p versus  z relation  in  the 
general  form  p = (k<£,  + kc/b)zn  is  inadequate  to  predict  the  pressure- 

sinkage  response.  Although  not  pointed  out  previously  in  this  paper,  it 
should  be  mentioned  that  our  analyses  to  date  are  based  only  on  the  experi- 
mental results  for  sinkage  in  the  range  0.  6n  to  2".  Perhaps  it  should 
be  added  that  Dr.  Bekker  has  not  claimed  that  his  equation  would  be  ade- 
quate in  the  entire  -£/b  region  we  have  studied. 

FUTURE  WORK.  While  much  has  been  learned,  there  is  clearly 
need  for  the  following: 

a.  Similar  laboratory  experiments  in  other  soil  media. 

b.  Experiments  with  greater  depth  differences  in  the  soil  bins  to 
assess  the  depth  effect,  if  any,  on  the  parameters  under  homo- 
geneous soil  conditions. 

c.  Revision  of  the  model  to  take  account  of  the  dimensions  of  the 
bearing  surface. 

d.  Improvement  of  the  model  to  cover  a ;wider  range  in  depth  of 
sinkage,  say,  from  at  least  0.5  to  5.0  inches. 
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TABLE  OF  MEANS  FOR  p*  VALUES 

(estimated  pressure  tor  one  inch  sinkage) 
summarized  by  plate  dimensions  (inches) 


SIZE  OF  PLATE 

P* 

COMBINED  MEANS  p* 

Circles 

2"  diameter 
4"  »' 

5.509 

8.222 

all  Circles 

6.865 

Rectangles 

1 x 4 

5.188 

Lengths 

(over  all 

Widths) 

1 x 6 

4.889 

4 

6.429 

l x 8 

4.701 

6 

6.378 

1 xIO 

4.737 

8 

6.296 

2x4 

2x6 

6.706 

6.671 

10 

6.141 

2x8 

6.420 

Widths 

(over  all 

Lengths) 

xIO 

6.1  70 

1 

4.928 

3x4 

7.392 

2 

6.492 

3x6 
3'X  8 
3x10 

7.576 

7.568 

7.515 

3 

7.5  12 
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Table  3 


ANALYSIS  OF 

VARIANCE  OF  p* 

VALl/ES 

(estimated 

pressure  for  one  inch 

sinkage) 

SOURCE  OF 
VARIATION 

DEGREES  OF 

freedom 

MEAN 

SQUARE 

Plates 

(13) 

8.4927 

Circles  vs  Rectangles 

1 

25279 

Between  Circles 

1 

22.0784 

Among  Rectangles 

(1  1) 

7.6181 

Widths 

2 

409546 

Lengths 

3 

0.4964 

Length  by  Width  6 

0.0668 

Plates  by  Depth 

(13) 

0.2290 

Experimental  error 

52 

0.1793 
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Table  4 


ANALYSIS  OF  VARIANCE  OF  n VALUES 

(estimated  from  six  replicates  grouped  to- 
gether for  fitting  the  pressure- sinkage 
equation  in  double -log  form) 


SOURCE  OF 
VARIATION 

Degrees  of 

freedom 

MEAN 

SQUARE 

Plates  (13) 

Circles  vs  Rectangles 

1 

0.0112 

Among  Circles 

1 

0.124  8 

Among  Rectangl-es  (it) 

Widths 

2 

0.  1 1 1 3 

Lengths 

3 

0.0228 

Lengths  by  Widths 

b 

0.0193 

Plates  by  Depth 

13 

00010 

Experimental  Error 


52 


00037 
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